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Signature of the rjNN configurations in coherent 7r° photoproduction on the deuteron 
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The photoproduction of a neutral pion on the deuteron is considered in the energy region around 
the rj threshold, where a bump-like structure was observed at very backward pion angles. Different 
dynamical aspects which may be responsible for this phenomenon are analysed within a theoretical 
frame which includes intermediate rjNN configurations. The results show in particular, that a 
three-body treatment of the rjNN interaction is of special importance. 
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I. INTRODUCTION 
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In various scattering and production reactions of pions on the deuteron, the cross sections exhibit a bump at 
■ backward angles for energies around the r] production threshold. In early discussions of the underlying dynamical 
mechanisms the main emphasis was put on the interpretation of these anomalies in terms of dibaryon resonances Q . 
. Later work, however, was mainly focused on a more natural explanation by relating this structure to the threshold 
' effect caused by the opening of the rj production channel. In particular, in a recent experiment on coherent tt° 
photoproduction such a bump was observed in the cross section in the backward direction around a lab photon energy 
t-H • Ej = 700 MeV 0. This nontrivial energy dependence was explained in in which the authors conclude that 
the mechanism, where first an rj meson is produced on one nucleon which then interacts with the second nucleon 
followed by 7r° production, is responsible for the enhancement around the rj threshold. In particular, according to 
the conclusion of y|, the appearance of the peak is nothing else than the signature of the 511(1535) resonance whose 
. excitation is believed to dominate rj photoproduction — > rjN. 

The purpose of the present paper is to perform a more detailed and extensive investigation of this phenomenon. 
. The main question we address here is the same as in |3j, namely what is the dynamics behind this enhancement. 
Although it might generally be clear that the observed structure is caused by the appearance of an rj meson in one 
way or another, the question about the underlying mechanism seems to be non-trivial. To be more specific, we would 
1 firstly like to note three main points which will be discussed separately, 
i— I \ (i) The wide peak in the jd — > n°d cross section can be a direct consequence of the cusp-like structure in the 
elementary amplitude — > tt°N in the Sn channel near the rj threshold. The case in point is the strong coupling 
between the irN and rjN states in the region of the 5n(1535) resonance resulting in a very pronounced cusp in 
the electric dipole amplitude E 0+ for pion photoproduction at E 1 w 710 MeV. The latter was also observed in the 
energy independent multipole analyses (see, e.g., energy independent solution in Ref. £j). Turning now to the process 
on the nucleon, which is bound in the deuteron, the elementary amplitude is expected to undergo an energy shift 
and a broadening of its structure due to the Fermi motion. Of crucial importance is the question to which extent 
this modification could affect the cusp in the elementary Eq + multipole and how prominent could be the Sii(1535) 
resonance in the reaction on the deuteron. 

(ii) The bump-like structure could also arise from an additional mechanism, appearing when the elementary photo- 
production process is embedded into the deuteron. Namely, one could expect that the anomalies arc, at least partially, 
caused by the three-particle unitary cut in the amplitude — > ir°d starting at the energy of the rj threshold. This 
cut arises because of the possibility to exchange a physical rj meson between the Six resonance excited on one of the 
nucleons and the second nucleon. The exchange mechanism is characterized by the pole which turns into the cut 
after the loop integration. Thus, the opening of a new physical channel leads to an additional contribution in the 
imaginary part of the amplitude reflecting a new inelasticity. Such a picture is typical for coupled channels and, if 
the corresponding dynamical equations are exactly solved, results in the three-body unitary relation. Because here 
we take into account only the leading term in the multiple scattering series, the whole amplitude does not fulfill the 
unitary relation. However, the r/NN three-body cut appears already at this level, and the question is, how does the 
amplitudes behave at the branch point. 

(iii) A sizable attraction in the rjNN system can lead to a strong correlation between all three particles. It is 
worth noting that in some work the rjNN interaction is predicted to generate even a bound state in the quasi- 
deuteron configuration (J W ;T) = (1~;0). Results provided by more sophisticated models [EEl show, however, that 
the fundamental rjN interaction is likely to be too weak for yielding binding of the r]NN system, so that only a virtual 
(antibound) state can be generated. In the context of the present discussion two points are relevant: (a) Although 
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the pole 'recedes' to the nonphysical region, it remains quite close to the zero energy. As a result, the virtual state 
can strongly influence physical processes involving an rj meson. Indeed, a variety of theoretical calculations and 
experimental analyses exhibit a strong rise of the rj production cross section just above threshold. This collective 
effect, in which all three particles participate, is naturally explained within a three-body model 6]. In this sense, the 
origin of the pole in the amplitude is not of crucial significance, and the existence of an rjNN bound state is not the 
necessary condition for an anomalous behavior of the rj production cross section. Important is only how far is the pole 
from the threshold energy, (b) Perturbative models, like the first order rescattering approximation noted in (ii), where 
only the leading order terms are kept in the multiple scattering series, are unable to reproduce the real dynamics of 
the ijNN system in the low energy region. The most simple explanation is that the corresponding Neumann series 
converges very slowly near the pole position jjj , so that the leading terms turn out to be a bad approximation to the 
whole series. From the last considerations we conclude that any realistic study of the role of rj mesons in the reaction 
jd — > TT°d should be based on a three-body approach to the r/NN system. 

Resuming now the qualitative discussion, we would like to note that all mentioned factors can come into play to 
form the observed characteristic bump in coherent tt° photoproduction close to the rj threshold. Important is the 
quantitative relation between the different mechanisms, which is the main object of the present paper. In the following 
we consider all three points separately with special emphasis with respect to their contribution to the resulting cross 
section. A brief description of the formal ingredients in Sect. [H] is followed by the discussion of our results for the 
differential cross section of jd — ► ir°d in Sect. Hill An appendix contains the listing of various formulas used for the 
calculation of the reaction amplitude in the impulse approximation. 

II. THE FORMALISM 

We start the brief description of the basic formal ingredients by presenting in Fig. ^ the diagrams which we 
consider in the theoretical analysis. The three combinations (a), (a)+(b), and (a)+(b)+(c) present three different, 
successively improved levels of approximation to the reaction amplitude, corresponding to the three points (i) through 
(iii) discussed above. We will refer to them as impulse approximation (IA), first order rescattering approximation, 
and three-body calculation, respectively. 




FIG. 1: Diagrams for the reaction 7c! — > ir°d included in the present work: (a) impulse approximation; (b) first order rescattering 
contribution; (c) additional contribution from three-body dynamics only in the s-wave (J lr ,T = CP, 1); (d) equation for the 
amplitude Ts 1± of photoproduction of the Sn resonance on the deuteron. 
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The starting point of our formalism is the impulse approximation which is completely determined by the elementar 
amplitude for 7iV — * 7r°iV and the deuteron wave function. As elementary amplitude we take the MAID analy; 
and for the deuteron wave function the parametrization of the Bonn-potential model (OBEPQ-version) [lOj with 
inclusion of the tensor component. The latter is certainly very important in the region of large momentum transfers. 
The general expression for the amplitude is given in the Appendix. 

For the additional contributions of the diagrams (b) and (c), we adopt the following simplifications. Firstly, in 
the pion exchange contribution to the second diagram only the Sn resonance was taken into account neglecting the 
contribution of other resonances. Although this neglect leads to an underestimation of the pion rescattering effect 
it should not strongly affect the quality of our results because of the following reasons: (i) since rescattering of 
intermediate pions is not related to the threshold effects it does not contribute into formation of the peak structure 
which we discuss here. What we can expect is only a smooth change of the cross section at backward angles; (ii) the 
role of pion exchange seems to be not very essential. In Ref. Q where this mechanism was calculated more precisely 
the corresponding effect in the second resonance region is less than 20 %. 

Secondly, as is indicated in Fig. ^ the three-body problem for the rjNN system was solved only for the s-wave state 
(J W ;T = 0"; 1). As is shown in Ref. @ it is the state of lowest orbital momentum which is mostly distorted by the 
multiple scatterings between particles. Other states with higher orbital momentum can be included perturbatively 
within the rescattering approximation (diagrams (a)+(b)). In order to avoid double counting we remove from the 
three-body amplitude those diagrams which possess the same topology as the ones already included in (a) and (b) 
shown in the same figure. 

For the elementary amplitudes appearing in other diagrams (b) and (c) of Fig. we assume that the photopro- 
duction of an rj meson as well as its interaction with nucleons proceeds exclusively via the excitation of the S , n(1535) 
resonance. As mentioned above, the same ansatz was adopted for pions. According to this assumption the t matrix 
of meson-nucleon scattering is given by the conventional isobar model 

t a N^/3N(p,p ;W) = — — — — — , a, e {tt, 77} , (l) 

W — M() — lj n — 2^ — Zjtttt 

where W denotes the invariant energy and S Q the various self energy contributions from the aN channels with 
a G {ir, 77, 7T7r}. The t matrix is determined by the bare resonance mass Mq and the parameters of the vertex functions 
g a , for which we take a simple Hulthen form 

9*(P)=9*U + jp) ■ (2) 
The contributions to the self energy from the various channels are expressed in terms of g a (p) as 

S (W) = — 1 9a{q)2 q2dq (3) 

a{ ' 2n 2 J W-E N (q)-u a (q)+ie2uj a { > 


for a G {7T, 77}. Since the double pion channel ttt:N is not explicitly included in the calculation, primarily because of 
its rather weak coupling to the 5*11(1535), we parametrize, following 11], the corresponding self energy in a simplified 
manner as a pure imaginary contribution proportional to the three-particle phase space 

i W-M N -2m w 

^Tvn = -7:!^ • (4) 

For the photoproduction amplitude we take the same ansatz as in (JTJ where one hadronic vertex function is replaced 
by the electromagnetic vertex g^jy for 7./V — * Sn which depends only on the invariant energy W and is parametrized 
in the form 

4 / v n 

V a n I — ) , for W > M N 



9-y P (M N + m T ) , else, 

g in (w) = -om 9lp {w). 



with q v — a/ (W 2 — (Mm + m^) 2 ){W 2 — (Mm — m„) 2 ]/2W. The isospin dependence of the 5n(1535) photo excita- 
tion amplitude was taken according to the relation |l2( 



cr(7P -> rjp) 
a(^n — » rjn) 



0.67. (6) 
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TABLE I: Parameters of the r/N scattering matrix in Eqs. Q through JSJ. The values of A a , 7^, and Mo are in MeV. 
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2.00 
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2.51 


404.5 4.3 


1598 


5.502 ■ 10" 1 


-1.923 ■ 10~ 2 


1.018 ■ 10" 1 


2.255 • 10~ 3 


-7.042 ■ 10" 3 



The parameters, appearing in the expressions Q through JSJ, are listed in Table [I] They are chosen in such a way 
that, on the one hand, the reactions — > aiV and 7r~p — > ryn are well reproduced in the Sn channel as presented 
in our previous work pH , Il4| . On the other hand, the chosen parameter set predicts for the rjN scattering length the 
value a V N — (0.5 + i0.3) fm which we consider as an approximate average of the various values provided by the rjN 
analyses. 

The method which we use to solve the three-body problem for the rjNN interaction is described in Ref. |6| and we 
refer the reader to this paper for the details. Here we only would like to mention that the key point of the method is the 
separable representation of the driving two-body interaction in the irN, rjN, and NN subsystems. The corresponding 
t matrix for meson- nucleon multiple channel scattering is given by the isobar formula For the NN sector, we use 
the separable representation of the Bonn potential as given in Ref. 0] for the 1 Sq and 3 Si configurations. 

As is well known, the separable ansatz makes it possible to reformulate the three-body problem in terms of two-body 
scattering between quasiparticles. As a consequence, the solution of the problem is given by an amplitude Tg^ of the 
effective transition 7c? — > NSn as presented in Fig.QJd). The needed physical amplitude — ► ir°d is then obtained 
through an additional loop integration (the first diagram in Fig. ^c)). 

III. RESULTS AND DISCUSSION 

In order to demonstrate the quality of the impulse approximation, we present in Fig. |2 the resulting pion angular 
distribution, and compare them with the available data from the compilation of . Although the agreement is quite 
satisfactory for the forward angles, the data are significantly underestimated in the backward region. We will return 
to this point at the end of the discussion. A comparison with the IA calculations of Ref. also exhibits differences 
which are, however, not very dramatic. For instance, our cross section at E~ f = 700 MeV shows quite a monotonic 
behavior and does not possess a deep minimum at very forward angles as obtained in Ref. It is intuitively clear 
that as long as small angles are considered where the two-nucleon effects are minimal (except for pion two-body 
absorption) the magnitude of the cross section should be mostly determined by the elementary amplitude. Therefore, 
the difference between the two theoretical results has to be ascribed primarily to the differences in the corresponding 
elementary operators, especially in the spin-flip part (see Eq. ljA4p ). which accounts for the dominant fraction of the 
forward cross section on the deuteron. 




B n [deg] B n [deg] 

FIG. 2: Angular distribution for jd — > 7r°d at two photon energies. The curves are the results of the impulse approximation 
(IA). The data are taken from the compilation in |16f . 

Now we turn to our main results on the additional interaction effects as presented in Fig. [21 The curves show the 
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FIG. 3: Differential cross section for -yd — > n°d as function of incident photon energy for two pion emission angles in the jd 
cm. system. The curves show the results of the impulse approximation (dashed), first order rescattering (dash-dotted) and 
three-body calculation (solid). The dotted curve shows the contribution of the rj rescattering term alone (diagram (b) in Fig.Ql 
where the resonance propagators were substituted by constants (see text). The corresponding values of the four- momentum 
transfer squared in fm -2 are shown on the top abscissa. 



predictions according to the different approximations discussed in the introduction. As one might expect, at forward 
angles the cross section shows very little influence of these effects. With increasing momentum transfer, corresponding 
to increasing pion emission angles, small internuclear distances come into play and thus corrections to the simple IA 
calculation from the two-nucleon mechanisms become more and more important. Furthermore, a bump in the energy 
dependence around E 1 = 650 MeV is clearly visible at very backward angles. Among other things, this fact can be 
considered as strong evidence that primarily the two-nucleon mechanisms are responsible for this phenomenon. 

This obvious statement does not, however, diminish the role of the single nucleon response. As one readily sees 
in the right panel of Fig. [21 some nontrivial structure, a shoulder, appears in the cross section near the r\ threshold 
already in the impulse approximation where the second nucleon is not actively involved. A more detailed analysis 
reveals, however, that the enhancement is not caused by the presence of the 5n(1535) resonance. Rather it is a 
combined effect of different terms in the MAID amplitude which we use here. 

On the other hand, we would like to note that the 5n(1535) resonance itself does in principle produce a slight 
shoulder close to the r\ threshold, which, as already mentioned in the introduction, is a signature of the cusp in the 
i?o+ multipole smeared out by the Fermi motion in the deuteron. However, this resonance contributes little to the 
coherent reaction on the deuteron, so that this cusp-like structure turns out to be invisible in the cross section. Hence, 
the slight enhancement observed in the otherwise monotonic behavior of the IA cross section is not related to the 
rjNN dynamics and should be ascribed to properties of the elementary pion production amplitude in this region. 

Now, when the first-order rescattering of the produced tt and r\ mesons is included, a peak structure clearly evolves 
as exhibited by the dash-dotted curve in Fig. [3] As was emphasized in the main origin of this effect is the presence 
of the 5*11 (1535) resonance in the diagram with r\ rescattering. Our results in general confirm this statement. It is, 
however, interesting that the new three-body unitary cut in the amplitude connected with the rjNN channel also 
make a slight contribution to the formation of the bump. Indeed, if the strong energy dependence of the amplitudes 
-fN -> T]N and r/N -> tt°N (see diagram (b) in Fig. ^) near the r\ threshold is 'neutralized' by keeping the 5*11 
propagator constant, the rescattering term alone exhibits slight enhancement (dotted line in Fig.[3J). 

Turning now back to the discussion of three-body effects in ir° photoproduction, the three-body calculation clearly 
predicts a much more prominent peak structure accompanied by a slight shift to lower energies. Thus the difference 
between the solid and the dash-dotted curves demonstrates convincingly the importance of the higher order terms in 
the multiple scattering series for the intermediate rjNN interaction. It is worth noting that in the 'ideal case' the 
measurement of the peak might be an indicator of the dynamical properties of the r]NN system. In particular, if a 
bound j]NN state would exist, it would appear in the irNN channel as a pure s-wave three-body resonance. The 
mass difference 2Mm + m v — Mr, where Mr is the mass of this hypothetic resonance, would then give the binding 
energy. If on the other hand the ijNN scattering amplitude possesses only a virtual pole, as is predicted by our model, 
the resonance peak should fall directly on the ij threshold. Nevertheless, it is obvious that the extraction of such 
information would require a quite thorough partial wave analysis (we need to consider only the s-wave contribution) 
and, moreover, a very precise energy resolution and a small statistical error, which could then allow one to fix the 
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FIG. 4: Comparison of our three-body results (solid line in Fig.|^ with CLAS preliminary data 0- 



resonance position. A similar scheme was already realized in the experiments aimed at a search for r\ mesic nuclei with 
3 He 0] an d heavier targets [l8| ■ But although the quality of the data permits one to make quite definite conclusion 
about the peak position, there are still problems with the interpretation of the measured angular distribution and the 
energy dependence of the corresponding cross section within existing theoretical models [l7| • 

A comparison of our results in Fig. 0] with preliminary data from Ref. |2( shows that the theory underestimates the 
observed cross section at backward angles by about one order of magnitude. The same discrepancy is noted for the 
results of Ref. 0] where the cross section is also far below the data in the same region. Yet, the authors of use 
an oversimplified operator for jN — > nN and it is therefore difficult to identify the prime reason of this drawback. 
Furthermore, similar underestimation may be exhibited in the work of Q, although the deviation is not as significant 
as in our case. 

The noted discrepancy might not be very surprising after all. We can expect that as long as the pions are produced 
on the individual nucleons and as long as the angular dependence of the elementary amplitude 7 AT — > ir°N is not 
varying strongly, the form of the differential cross section is mainly governed by the deuteron form factor. In order 
to locate the appropriate portion of the form factor which enters in the kinematical conditions, we present on the top 
x-axis in Fig. [J] the corresponding four momentum transfer squared. As one can see, its characteristic values range 
between 15 and 35 fm -2 , where the form factor as seen in electron scattering exhibits a sizeable sensitivity to higher 
order mechanisms like tt and p MEC's (see, e.g., Ref. [lflj)- 

Therefore it is very probable that additional two-nucleon mechanisms, not included in the present calculation, will 
make sizeable contributions to the backward cross section. For example, meson exchange currents were found to be 
quite significant for n + photoproduction on 3 He in Ref. [23] • In this case, the effect is associated with double meson 
photoproduction followed by absorption of one of the mesons on the spectator nucleon. However, if such mechanisms 
are really quite important, their inclusion would tend to diminish the effect discussed above, i.e. the relative size of the 
peak as its manifestation. In this case, we will have to face again the same question about the origin of the structure 
observed in the ir° photoproduction. 



IV. CONCLUSION 



We have discussed various aspects of the influence of rj degrees of freedom on coherent ir° photoproduction on the 
deuteron. Two interrelated features determine the significance of the rj meson in this reaction. Firstly, the strong 
coupling between ttN and rjN states in the energy region of the S , n(1535) resonance leads to a significant admixture 
of the i]NN configuration to the NSn intermediate states. What is more important, in contrast to pions, slow 77 
mesons interact strongly with nearby nucleons. As predicted in a variety of investigations, this dynamical feature 
leads to highly correlated r/NN states, which should manifest themselves in the s-wave part of the ir° photoproduction 
amplitude. 

Both factors lead to the appearance of a pronounced bump in the energy dependence of the backward differential 
cross section around the 77 production threshold. This effect was analysed in Ref. Q on the basis of a theoretical 
model where in addition to the mere impulse approximation 7r and 7/ rescatterings were included. In contrast to the 
conclusion of [3j, we find that already in the impulse approximation a shoulder appears in the cross section around 
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-E 7 = 700 MeV. However, as is explained in the present work, this effect does not have a deep physical significance, 
and is likely to reflect the special structure of the pion production amplitude used here. 

The inclusion of first-order rescattering and, finally, all terms in the multiple scattering series within a three- 
body model shifts the peak position and make it significantly more pronounced. In particular, it was shown that a 
three-body treatment of the tjNN interaction is of special importance for the understanding of the reaction dynamics. 

In general, according to our results, the physics behind the bump structure in the cross section for 7c? — ► ir°d may 
be much more complicated than was presented in Ref. 0. Among other things, we are very sceptical about the 
possibility to extract a model independent information on the fundamental rjN interaction from the (7, tt°) reactions. 

The present calculation as well as those presented in Ref. Q can be considered as a natural explanation of the 
experimental results reported in 0. However, the status of these findings is quite unclear, because the theory strongly 
underestimates the data in the relevant angular region, making any quantitative analysis of the observed cross section 
impossible. Further theoretical investigations in this field are certainly needed. 
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APPENDIX A: IMPULSE APPROXIMATION FOR THE AMPLITUDE jd -> iv°d 

Here we give a brief outline of the impulse approximation of the T-matrix of coherent pion photoproduction on the 
deuteron in the cm. frame 

7 (w 7 , k, A) + d{E d , -k) -> ttK, q) + d{E' d , -q) , (Al) 

where energy and momenta of the participating particles are given in the parentheses, and A stands for the circular 
photon polarization index. In the impulse approximation, the amplitude T mm i\ for the transition between the target 
states with spin projections m and ml on the z-axis, chosen along the photon momentum, 

44s tin'iP 1 ) t^(k,p h qj f )4> m {p) ■ (A2) 

with standing for the corresponding elementary amplitude 7 A — > n°N. Furthermore, the vectors pi and pf denote 
initial and final momenta of the active nucleon in the deuteron, for which we have pi — p— k/2 and Pf = p '— q + k/2, 
and <p' = p + (k — <f)/2 denotes the relative momentum in the final deuteron state. 
For the deuteron wave function we use the familiar ansatz 

<t>m{p)= X] H (Lm L lm s \lm)u L (p)Y LmL (p)x ms (o , (A3) 

L=0,2 m L ms 

where the last two terms denote spin and isospin wave functions, respectively. 

For nuclear application it is convenient to split the amplitude into spin independent and spin-flip parts (the 
index A is omitted for convenience in the following expressions) 

t wl = K + iL ■ a . (A4) 
Then, using standard angular momentum algebra, the reaction amplitude i|A2fl can be put into the following form 



Tw = A^^(-l) M AV2ATl(lmA-M A |lmO £ / 7^3 [{f I \ } O^V) ® Y lL] (P)}[1 A K 

A=0 LL'=0,2 \ n ) { J 

_ 3 (All) 

- (-1)W\^TT 1 1 1 \{{Y^'\p')®yM(P)}W®lM} [ $ u L ,(p')u L (p) . (A5) 

i=o [ILL') J 

As for the isospin structure, it easy to understand that from all three amplitudes in the isospin decomposition of 
the elementary operator for pion photoproduction with Cartesian index a=l, 2, 3 pl| 

U 7 = M^T a + M(-^[r a ,T 3 ] + M {+) 5 a3 , (A6) 
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only can contribute to the coherent process on the deuteron. 

Using standard normalization of particle states, the cross section related to the same cm. frame reads 



da q E d E' 



dn w 7 UirW) 2 6 

7 v ' mm'X 



\Y,\ T ^x\ 2 - (A7) 
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